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Abstract. We prove that the diagonal of the transition probabilities 
for the c?-dimensionaI Bessel processes on (0, 1], reflected at 1, which we 
denote by p^(t,r, r), is an increasing function of r for d > 2 and that 
this is false for d =2. 



1. Introduction 

The following conjecture of Richard Laugesen and Carlo Morpurgo arose, 
as communicated to us by R. Laugesen, in connection with their work in 
[12] on conformal extremals of zeta functions of eigenvalues under Neumann 
boundary conditions. While this may be the first time the conjecture ap- 
pears in print, the problem seems to be well-known. 

Conjecture 1.1. Let IB he the unit hall in IR^, d>2, and let p^{t, x,y) be 
the heat kernel for the Laplacian in IB with Neumann boundary conditions. 
Equivalently, p^{t,x,y) gives the transition probabilities for the Brownian 
motion inlBwith normal reflection on the boundary. Fixt > 0. The (radial) 
function p^{t,x,x) increases as \x\ increases to 1. That is, for all t>0, 

(1.1) p^{t, Xi,Xi) < p^{t, X2,X2), 

whenever < < \x2\ < 1. 

Of course, the same conjecture makes sense for d = 1. For this, see 
Remark 5.4 in §5. 

We should observe here that for the Dirichlet heat kernel in IB, the op- 
posite inequality is true. That is, the diagonal of the Dirichlet heat kernel 
decreases as the point moves toward the boundary (see Proposition 5.2 in 
§5 below). 

The conjecture is closely related to the hot-spots Conjecture of Jeff Ranch 
which asserts that the maxima and minima of any eigenfunction ipi corre- 
sponding to the smallest positive Neumann eigenvalue /.ti of a convex planar 
domain are attained at the boundary, and only at the boundary, of the do- 
main. Indeed, if we denote the volume of the unit ball IB in IR' by Ud, the 
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eigenfunction expansion of the heat kernels gives that 
(1.2) pN^t,x,x)^- + e-^-'\^i{x)\\ 

and this is uniform in x for t large (see [15]). We refer the reader to [2], 
[8], [4], and references therein, for more on the hot-spots conjecture and 
for the use of heat kernel expansions and transition probabilities for that 
problem. Of course, for the unit ball the hot-spots conjecture follows easily 
from the explicit expression of (pi as a Bessel function. However, a more 
general Laugesen-Morpurgo Conjecture can be stated where the connection 
to the hot-spots conjecture is more meaningful, see Conjecture 5.1 below. 
Surprisingly the hot-spots conjecture is open even for an arbitrary triangle in 
the plane. But perhaps even more surprising is the fact that the Laugesen- 
Morpurgo conjecture is open even for the unit disk in the plane. 

The Neumann heat kernel p^{t,x,y) gives the transition probabilities 
for the Brownian motion reflected on the boundary of the ball and hence 
the use of probability for this problem (just as in the case of the hot-spots 
conjecture) is very natural. The Brownian motion in the ball has a skew 
symmetric decomposition in terms of a Bessel processes (the radial part) and 
spherical Brownian motion running with a clock that depends on the Bessel 
processes (see for example [<i]). That is, let be reflected d-dimensional 
Brownian motion (RBM) in the ball IB and let Rt be the d-dimensional 
Bessel process in the interval / = (0, 1] reflected at 1. Then B4 is the 
radial part of W^. That is, Rt = \WP\. Let pf{t,r,p) be the transition 
probabilities for Rt in the interval / = (0, 1]. We will often refer to this as 
the heat kernel for | W^\ . The main result of this paper is the following 

Theorem 1.2. Suppose d > 2. Fix t > 0. The function pf(t,r,r) is 
increasing in r. That is, 

(L3) pf{t,n,ri) <pf{t,r2,r2) 

for < ri < r2 < 1. This monotonicity property fails if d = 2. 

It is known (see [14], page 415) that the transition probabilities (heat 
kernel) q{t, r, p) for the free 2-dimensional Bessel process is given by 

(1.4) Q(t,r,p)=t-V-T/o(!^), 

where Iq is the modified Bessel function of order 0. The function q[t, r, r) 
is not increasing. In fact, as t — > this function has a "tall bump" moving 
toward 0. Since the reflected process is equal to the free process before the 
first reflection, it seems reasonable to expect that pfit, r, r) is not increasing 
for small values of t when d = 2. A rigorous proof of this fact will be given 
below. On the other hand, the function q{t, r, r) is non-decreasing for d- 
dimensional Bessel processes when d >3 and therefore one may expect the 
monotonicity property to hold for pf{t,r,r) for such d and Theorem 1.2 
shows that this is indeed the case. 
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Our strategy in this paper is to replace the reflected Bessel process by a 
random walk and obtain the result for this walk. We then show that under 
the appropriate scale the random walk converges to the reflected Bessel 
process. 

The paper is organized as follows. In §2 we introduce the random walks 
and prove the analogue of Theorem 1.2 for these. §3 contains the proof of 
the convergence of these random walks to their continuous counterparts and 
§4 gives the proof of theorem 1.2. Finally, the last section contains some 
conjectures related to our result and illuminates the connection between 
the Ranch hot-spots Conjecture and the Laugesen-Morpurgo Conjecture 
further. 



2. The random walk 

In this section we introduce the random walk which we will use later to 
approximate the reflected Bessel process. We will use the following notation. 
For any x S W^, \x\ denote the length of the vector x. For d> 3 and x 0, 
set 



and 



TTt \ 1^1 + 1 

U (X) = — j — j — X 



D{x) = ^ x. 
\x\ 



Consider the two sets 

C{x) = |y G M"^ : \y\ = + 1 and y - D(x)_Lx| , 

and 

5 = |y G R'' : |y| G N and |y| > ^ - l| . 

Note that C{x) is a (d — 2) -dimensional sphere with center at D{x) and 
orthogonal to x. We now define two random walks as follows. 

Definition 2.1. Let be a random walk on S with the following transition 
probabilities 

(1) p{l,x,U{x)) = \, 

(2) p(l,:E,D(x)) = i-g,/or|x|>^^, 

(3) p(l,j;, A) = ^^/^^.(A), where A C C{x) and /i^ a uniform probabil- 
ity measure on C{x), for \x\ > 

(4) p{l,x,A) = ^iJx{A), where A C C{x) and fix a uniform probability 



measure on C{x), for even d and \x\ 



d-2 

2 • 
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Observe that if d = 2A; + 1, then > k and p{l, x, D[x)) = if |x| = k. 
li d = 2k, then \x\ > k — I and x, D{x)) = ^ if = k. Hence, we need 
the additional points |x| = A; — 1 as in the last case in the definition. 

We will also be concerned with the radial part of the above random walk. 

Definition 2.2. Let Yn = Then Yn has the following transition prob- 

abilities 

(1) p{l,m,m + l) = ^ + ^, form>^, 

(2) p{l, m,m + 1) = 1, for even d and m = ^ — 1, 

(3) p{l, m,m-l) = \-^, form>^. 

The above are the "free" random walks. For our purpose, we define the 
reflected versions of these walks. 

Definition 2.3. Let be the random walk on 5'n{|x| < with transition 
probabilities pN^'-.x^y) = p{l, x,y) for \x\ < N and 

(1) pn{1, x,x) = ^ + ^ for \x\ = N, 

(2) pn{1, X, D{x)) = \-^ for\x\=N 

and denote by Y^ = \X^ \ its radial part. 

We will use and Y^ to approximate and Rt = \W^\, respec- 
tively. We will now show that p]y{n,m,m), the transition probabilities for 
the random walk Y^ , are increasing with m for any fixed n. The following 
property will be useful for this purpose. 

Definition 2.4. We say, that a random walk has the nondecreasing loop 
property if for any m its transition probabilities have the property that 

p{\, m,m + m + 1, m) 

is nondecreasing with m. In case of a reflected walk we also require that 
pj\f{l, N, N)'^ (the loop at the reflection point) be larger than or equal to 

p{l, m,m + l)p(l, m + 1, m), 

for every m < N . 

Lemma 2.5. The random walks Y^ and Y^ have the nondecreasing loop 
property. 
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Proof. In the case of Yn we have for m > {d — l)/2 

/I d-l\ /I d-1 \ 
p{l,m,m + lMl,m + hm)=^- + ^ji^- " ^^^^J 



f2.ll 



_ (2m + d- l)(2m + 2-d+l) 

lQm{m + 1) 
_ 4(m^ + m) - (d- 3) 

16(m2 + m) 
_ 1 _ (d- l)(d-3) 
4 16(m2 + m) 

Thus the left hand side is nondecreasing for m > When d = 3 this 

quantity is constant. 

Next we take m = | — 1 for d even. In this case 

, , / d d\ ( d d \ 1 

(2-2) p 1,--1,-U =—. 



2 '2/'^V'2'2 / 2d 



From the general case 



,dd\(dd\ 1 {d-l){d-2,) 

p l,-,- + lU 1,- + !,-^ - 



2'2 7'^V'2 '27 4 4d2 + 8(i 

6(i-3 



> 



4d2 + 8(i 
2d + 4 1 



Ad^+M 2d' 

since d > 4. This completes the proof that Yn has the non-decreasing loop 
property. 

In the case of Y^ we are left with reflection point loop, i.e. p(l, N, N)"^. 
But this is larger then 1/4, hence we have a nondecreasing loop property for 

yN □ 

Proposition 2.6. Fix n. Then pi\i (12,171,171) is increasing in in. 

Proof. To prove this we fix n and consider each possible path from m to m in 
n steps. The proof will be completed if for each of them we can find a unique 
path from m + ltom + linn steps that has a larger probability. Towards 
this end, let Pm = {"m = h,l2,--- Jn-i,ln = rn} be a path for 1"^. Let 
ki = infjfc : = N} and ^2 = supjA; : l/. = N}. If this path never touches 
then we can take the path Pm+i = {m+1 = /i + 1, Z2 + I, • " " 1 + 1 = m+l}. 
Since both paths start and end at the same point, they are (possibly after 
rearranging) a sequences of loops. By the nondecreasing loop property for 
Y^ proved in Lemma 2.5, the probability for the path Pm+i is larger. 

Up to now we have only used those paths starting at m + 1 that do not 
remain at N. That is, those paths which move to — 1 immediately after 
hitting N. This is true, since all the paths Pm we considered above never 
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touched N, their maximum can be iV — 1 and the walk has to move to iV — 2 
from there. 

Suppose that Pm hits N at time ki < N. Then ^2 < ^- Let Pm+i = 
{m+1 = li + 1,12 + 1, •• • ,lki-i + l,lki, - ■ ■ ,^fc2'^fc2+i + l) • • • ,in + l = m+1}. 
The idea is to shift the parts of the path before and after hitting for the 
first and last time. We have to show that such correspondence of the paths 
is one-to-one. 

We have to look at the parts of Pm before, after and in between the 
hitting times, separately. First notice that for Pm+i the step from A;i — 1 to 
ki is the first time the path remains at N. Similarly k2 to A;2 + 1 is the last 
time the path remains in N. Hence if two different paths Pm have different 
times ki or k2, then the shifted paths Pm+i will also be different. Hence the 
only possible paths with the same corresponding shifted paths must have 
the same hitting times ki and k2- 

Note that it is important in this proof that the walk cannot stay at any 
point other then the reflection point A^. If we allow p{l,N — 1, N — 1) to 
be non-zero, then the path Pm+i obtained from the path Pm that never 
touches may remain at the point N even before the time ki. This would 
invalidate the above reasoning. 

If two paths Pm are different at any shifted point (before ki or after k2, 
then the same is true for the corresponding paths Pm+i- Finally, if both ki 
and ^2 are the same for two different paths Pm and they are the same before 
ki and after A:2, then there must be a difference between ki and k2- But this 
part of those paths is not changed in Pm+i- Therefore the correspondence 
between Pm and Pm+i is one-to-one. 

The last thing to check is that the probability of the corresponding Pm+i 
path is larger. Since the part between ki and k2 is exactly the same, we can 
disregard it. What is left is just a sequence of loops (after rearrangement). 
Hence by the nondecreasing loop property (Lemma 2.5), this completes the 
proof. □ 

As a corollary to the above argument we get 

Corollary 2.7. Fix n. For any m' < m, p > and m + p < N we have 
p]\f{n,m,m') < p]y{n,m + p,m' + p). 

The proof is almost identical. The only difference is that each path can 
be decomposed into loops plus additional transitions from m to m' (all of 
them toward 0). But since pN{^,Ti,n — 1) = | — is increasing with n, 
the path shifted by p will have larger probability. □ 

3. Convergence 

Proposition 3.1. The sequence {-^Xj^a^j} converges weakly to the reflected 
Brownian motion as iV — > oo. 

The proof of this fact is essentially the same as the convergence proof 
in [10]. First we need to establish the existence of a weak limit of the 
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process that interpolates ;^Xj^2j] linearly. That is, of the continuous 

process that equals j^X^2£^ at the times of the jumps of the process and is 

linear in between. The process Z^^ converges weakly by Holder continuity 
and Prohorov's theorem (see [10].) Our main goal here is to identify this 
weak limit as REM on IR To accomplish this we use the submartingale 
characterization of the reflected Brownian motion (see introduction in [Hi]). 
More precisely, the RBM in IB is the only stochastic process starting from 
x G IB such that for any / € C| (IB) with positive normal derivative at each 
point of the boundary of IB, the process 

(3.1) fiW^)- f M{Wf)ds 

Jo 

is a submartingale. Hence, to prove that Z[^ weakly, it is enough to 

show that 

(3.2) hm mf E (/(Zf ) - /(Zf ) - ^ £ A/(Z„^)dn) > 0, 

for all such functions /. Here and in the sequel, A denotes the Laplacian in 
IR^. 

First we will calculate an expectation of the single jump of ;^^[^2j] ■ Note 
that by the definition, Z^ = jjX^2f-^ at the jump times. 

Lemma 3.2. Let Un = n/N'^ be the points where the process jjX^2t] makes 
its jumps. Then 

7N \ f(r7N\\_-Eyx( 1 Af^r^N' 



[f{ZZ,J - f{ZZ) ) = (^A/(Z- ) + o(iV-) 
(3.3) +0(iV-^)l{|z-|=|#} 

where dif denotes the outer normal derivative of f on {y : \y\ = \x\}. 

Proof. Let A{x) = E^^(/(^Xf ) - /(x)). Note that if the starting point 
for the process Z^^ is x, then the corresponding starting point of is Nx. 
By the strong Markov property for X'^ and the definition of Z^ , 

E^ - /«)) = E^^- (/ - / 



(3.4) = E^^E^ 



JV 



Let Hx be the uniform probability measure on C{Nx)/N. For 1 > |x| > 
\x\ = k/N (the states of the rescaled process) and for any function 
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/ e cl 

(3.5) 



Let di denotes the outer normal derivative to the sphere {y : \y\ = \x\} 
at the point x. Let also d^i = didi. We have 

(3.6) 

Mx) = \ [^d,f{x) + ^9?i/(x')) 

+ ^) (-^d,f{x) + -^d!j{x") 



2 4N\x\J V 2iV2 
4A^kl Jc{Nx)/N 2 

+ JnTi I [(^1 - + ky - ^) • V]V(^(y))^/^x(y), 

4A^kl Jc{Nx)/N 2 

since C(x) is a sphere on d — 1 dimensional hyperplane orthogonal to x and 
centered in D[x). We also have yi — xi = —1/N hence the first order terms 
cancel and we get 

(3.7) . 

+ / [{y-x).V]'f{z{y))df,Ay). 

8N\x\ Jc{Nx)/N 

For X G 5(0,1), the function / has uniformly continuous second order 
derivatives. Hence the error in the Taylor expansion is uniformly bounded. 
If we denote the derivatives in the directions tangent to the sphere {y : \y\ = 
\x\} at X by di, i > 2, and di as before denotes the outer normal derivative 
to the sphere {y : \y\ = \x\}, then 

/ [{y-x)-Vff{z{y))dMy) 

JC{Nx)/N 



(3 8) = / [{y-x)■VYf{x) + o{\y-x\')d^l^{y) 

^ Jc[Nx)/N 



n 

' 'X\ 



y](2;i - yifdlf{x)dfia:{y) + o 

C(Nx)/N .^^ 



N 



since \y — = 4(|x| + 1/N)/N and mixed derivatives disappear due to the 
symmetry of C{Nx)/N. 
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If i = 1 in the above sum, then (xi — yi)^ = For each i > 2 the 

integral has the same value due to the symmetry of C{Nx)/N . Also 

n /III 

(3.9) 



Hence 



i=2 



A{x) = + o(A-2) - + o{N-M 



(3.10) +'^~^ 



m\x 
1 



A/(i) + o(JV-^) 



2 

2 



2A2 

since 2A|a;| > d — 1. 

Now we have to consider two special cases. First suppose d is even and 
\x\ = ^0-. Then 

(3.11) Aix) = 0(A-2) = ^Af{x) + 0{N'^). 

Note that the error is uniform, just like in the first case. 

The remaining case is the reflection circle. That is, the points \x\ = 1. 
These points correspond to the times when |Xj^2j]| = A. That is, when 

the walk is on the boundary. Hence we have to use the transition steps 
from Definition 2.3. Now, 

(3.12) 

= -CA,a,/w + ^A/w + 0{N-^). 

As before, di is the outer normal derivative to the sphere {y : \x\ = \y\}. 
Hence —di becomes a normal derivative along the reflection direction if x is 
on the boundary (|x| = 1). By the definition of / the first order term above 
is positive. 

Combining all the cases we obtain 

(3.13) = 2^^^(^) + + 

+ (-c^ai/(x)+0(A-2))l{|,|^i|. 

This and (3.4) give the assertion of the lemma. 
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□ 

We are now ready to prove (3.2). Summing the expressions from the 
above lemma over s < u„ < t yields 

(/(Zf ) - /(Zf )) = (^1^ ^A/(Z„^)dn) 

'ct,{s,t))0{N~^) 
y 2N / 



(3.14) 



+ E^ 



S<Un<t 

+ W (£f(s,t))0(iV-2) + o(l), 

where {s,t) denotes the local time (number of visits) of Z'^ on the set 
l^l = a between times s and t. Note that the term involving —dif is always 
positive, since by the definition / has a positive derivative in the reflection 
direction. In order to finish the proof of (3.2), we need the following lemma 

Lemma 3.3. Let Ly^ = ^ |n : = y, n < N^^ be the local time at y. 
Then 

(3.15) Ey{L^') = o{N'). 

Moreover, ^^{Ly'^) is increasing with y. 

We need to reduce the local time of the process to the local time of 
the process . Since t is fixed, the expectation of the local time {s, t) is 
comparable (with constant depending on t but not on N) to the expectation 
of same on the interval to 1. More precisely, by the strong Markov property 
for any a < 1 we have 



W{L^ {s, t)) < E"(£^ (0, t)) =E^{C^ (0, 1) + E^i {C^ (0, t - 1)) 

(3.16) 



<E" (£^(0,l)+E-(£^(0,t-l))) 



E"(£^(0,l))+E"(£^(0,t-1)) 
<•••< rtlE"(£^(0,l)), 
where [t] is the smallest integer bigger or equal to t. 



But, the local time (0, 1) of Z^ is the same as the local time Lf 



of Y^ . Hence, by Lemma 3.3 both error terms involving the local time 
in (3.14) are negligible. This ends the proof of (3.2). Hence the process 
-^Xj'J^2^] converges weakly to . It follows, that jqY^2£i^ converges weakly 

to Rt- The monotonicity of p^{t,x,x) will follow from the monotonicity of 
the approximating random walk, as we shall show in §4 below. 

Proof the Lemma 3.3. The idea of the proof is based on the following well 
known facts for the random walk Z„ on Z, Zq = 0, P{Zn+i = Z„ + 1) = 
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P{Zn+i = Zn — 1) = 1/2. It is a standard fact (which follows from the 
reflection principle) that the number of paths of length 2n satisfying 

{Zi > 0, Z2 > 0, . . . , Z2„_l > 0, Z2n = 0} 

is 

1 /2n - 2\ 
n\n — I J 

Hence, 

P{Zi > 0, Z2 > 0, . . . , Z2„_l > 0, Z2n = 0} = - 2-2". 

n\n — 1 J 

Next we show that E'^ilf) = o{N^). Of course, {hf) > E''{Lf) 
so it is sufficient to show that E^ {L^) = o{N'^). Consider the sequence of 
stopping times Rq = 0, Rk+i = inf{m > : (m) = N}. We have 

E^{Lf) = E^{max{k e N : Rk < N^}). 
We have Rk = X]j=i(-^i ~ Rj-i)- {Rj — is a sequence of i.i.d. 

random variables with Rj — Rj^i = -Ri. Let {Si}°^^ be a sequence of 

i.i.d. random variables with Si = Ri = inf{m > : {m) = and 

Ti = S,^2[N/A]. We have 

E^{Lf) = E^{max{k e N : Si + . . . Sk < N^}) 
(3.17) < ^^(max{A; G N : Ti + . . . < N^}). 

Note also that N - 2[A^/4] > N/2 > 2[N/4\. We may and do assume that 
is large enough so that 2[N/4] > {d-l)/2. 

Our next aim is to estimate E^T\. Let 2n < 2[A^/4]. The number of 
paths of length 2n satisfying 

{Fo^ = N, Yf <N,Yi' <N,..., Yl^, < N, Y^ = N] 

is 

1 f2n - 2\ 
n\n — 1 J 

Using Definitions 2.2 and 2.3 and the fact that m> N -2n> N - 2[Af/4] > 
2[7V/4], we obtain 

P^(5i = 2n) = P'^iY,^ = N,Yl' <N,...,yZ_^<N,yZ = N} 
- n\n-l) \2 ~ 8[N/4]) 

Observe that 

V2 8[N/A]) - 22" V 4[iV/4]y -2'^"-' 

We now adopt the convention that c = c{d) > is a positive constant which 
can change its value from line to line. 
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By Stirling formula 
1 /2n-2\ 1 (2n-2)! 



n V n - 1 y 22" n((n - 1)!)222" 

(2n - 2)2«-2e-{2n-2)^27r(2n-2) 



> 



> 



n{n - l)2("-i)e-2{«-i)27r(n - l)e2/(i2(n-i))22r, 
c 



^3/2 

Hence P^(5i = 2n) > c/n3/2. j^. foUows that 

E^Ti > 2nP^{Si = 2n) > cN^'^ . 

n=l 

Now we will estimate (3.17). Note that Ti > 1 so max{A; EN: Ti + 
■ ■■Tk< iV^} < N"^. 

Let M = [2N'^/E^T{\ + 1 so that ME^Ti - N'^ > N^. Note that 
M < cN^/"^. We have 



Ti + ...Tk<N^}) 



Ti + . . .Tk} < N^; Ti + . . .Tm > N"") 
Ti + ...Tk} <N^; Ti + ...Tm<N^) 



= E^{max{k G N 

+ ^^(max{A; G N 
(3.18) < M + N'^P'^{Ti + ...Tm <N'^). 

We have 

P^{Ti + ...Tm < N^) 
< P^{\Ti + ...Tm - ME^Ti I > ME^Ti - iV^) 
E^{\Ti + ... Tm - ME^Til'') ^ ME^Tf 
^ ' - {ME^Ti-N^f - ■ 

Recall that M < cN^I'^ and Ti = Si A 2[N/^] < N. It follows that E^T^ < 
A^2 and the last expression in (3.19) is bounded from above by 

< 



m - AfV2- 

By (3.18) we obtain E^{L^^) < cN^'^ which gives E^{Lf) = o{N^). 
The monotonicity for the local times follows from 

E,(Lf ) = E,(J] 1{Y,^=.}) = Y.^.iYi' = x) = Y,PN{k,x,x) 

k=Q k=0 k=0 

k=0 

where the inequality follows from the heat kernel monotonicity obtained in 
Section 2. □ 
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4. Proof of Theorem 1.2 

First we give the proof that if d = 2 then p^{t, r, r) is not increasing for 
smah enough times t. 

Let P C be a convex polygon and denote its Neumann heat kernel by 
Pp{t,x,y). It is proved in [9] that 

(4.1) lim4^ = l 

p{t,x,y) 

uniformly in x,y £ P, where p denotes the heat kernel of the free Brownian 
motion in R^. In addition, [9], also proves that if Di is a convex domain 
whose closure, Di, is contained in the convex domain D2, then there exists 
a to sufficiently small such that 

(4.2) PD,it,x,y)<PDS^,x,y), 

for all x,y £ Di and < t < to, where to depends only on the distance 
between dDi and dD2. By taking two polygons Pi and P2 such that Pi C 
IBc IBc P2 and combining (4.1) with (4.2) we see that 

(4.3) lim4^^ = l 

*-o p{t,x,y) 

uniformly in \x\ < ^ and \y\ < ^. 
Since for any x £JB , 

(4.4) P^-{|I^,^| G (a, 6)} = r r pp^{t,x,pe^')de 

Ja Jo 

and the reflected Bessel process is the radial part of the reflected Brownian 
motion, we have 

(4.5) pf{t,r,r)= rp^{t,r,re'^)de. 

Jo 

and similarly, 

/•2tt 

(4.6) q{t,r,r) = / rp{t,r,re''^)d6. 



Jo 

Let e > and r,p < 1/2. Using (4.3) we can pick t{e) such that for 
t < t{e) 

(4.7) (1 - £)p{t, r, re'^) < p^{t, r, re^^) < (1 + e)p(t, r, re^^) 

(4.8) (1 - e)p{t, p, pe'') < pg{t, p, pe'') < (1 + e)p{t, p, pe''). 
From this and the integral formulas for the Bessel heat kernels above 

(4.9) (1 - e)g(t, r, r) < pf{t, r, r) < (1 + e)q{t, r, r) 

(4.10) (1 - e)qit, p, p) < pfit, p, p)<il+ e)qit, p, p). 
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From (1.4) we have 

(4.11) q{t, rVi, rVi) = t-^/^e-''' Io{r^) 

Set ^Q{r) = re~^^Io(r^). Using tables of the Bessel functions one can check 
that $o(l) ~ 0.4657 and that $o(2) ~ 0.4140. Hence $o is not nondecreas- 
ing. Let r < p be such that $o('") > ^o(p)- This is the same as 

q{t,ry/i,ry/i) > q{t, pVi, pVi). 

Pick £ small enough to have 

(4.12) (l-e)$o(r) > (l + e)$o(p). 
Now for any t we have 

(4.13) (1 - e)q{t, rVt, rVt) > (1 + e)q{t, pVt, pVi). 

Set ri = r\/t and r2 = pVt so that ri < r2. If we take t small enough to 
have t < t(e) and r2 < ^, it follows from (4.9) and (4.10) that 

(4.14) pfit,n,ri) > {l-e)q{t,ri,n) > (1 + e)g(t, ra, rg) > (t, ra, rg). 

This completes the proof of Theorem 1.2 when d = 2. 

Now we turn to the case d > 3. Fix x with r = |x| < 1. Let fe{x) = 
X[r-e,r] , £ > 0. We have 

(4.15) E-' (f, (lY,^.^) = (Iy.^.. e [r - e,r] 



The event above consists of transitions from r to some points to the left of 
r. Hence by Corollary 2.7 this probability is increasing in r. By the weak 

convergence of (^-^^^21^ to {R^), we have 

(4.16) (a (^^[^2,])) - E^(/,(i?f )) = G [r - e,r]). 

Since the limit of increasing functions is nondecr easing, we have that for 
arbitrary £, 

(4.17) / pf{t,ri,p)dp< pf{t,r2,p)dp,iiri<r2. 

Jri—e Jr2~e 

Suppose that pf {t,ri,ri) > pf{t,r2,r2) for some ri < r2. By the conti- 
nuity of the heat kernel there exists e > such that 

(4.18) inf pf{t,ri,p)> sup pf{t,r2,p) 

pe[ri-£,ri] pg[j.2_e,r-2] 

Hence, 

(4.19) pfit,n,p)dp> r pf{t,r2,p)dp. 

Jri~e Jr2~e 

But this contradicts (4.17). Thus pf{t,r,r) is nondecreasing in r. 
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For any t > 0, the function pj[t,r,r) is a real analytic function of r G 
[e, 1], for any e > 0, since it is the diagonal of the heat kernel of an operator 
with real analytic coefficients. Thus, if it is nondecreasing then it must be 
strictly increasing. This completes the proof of Theorem 1.2. □ 

5. Further remarks 

As mentioned above, the Laugesen-Morpurgo conjecture implies Rauch's 
hot-spots conjecture for the disk. Of course, as already also mentioned this 
is a trivial observation since for the disk the Neumann eigenfunctions are all 
explicitly known (Bessel functions) and the hot-spots conjecture is trivial by 
"inspection". This observation, however, leads to a more general problem 
for planar convex domains where the connection to the hot-spots conjecture 
is more meaningful. 

Conjecture 5.1. Suppose 0, is a bounded convex domain in the plane which 
is symmetric with respect to the x-axis. Letp^{t, z, w) be the Neumann heat 
kernel for Q. Thenp^{t, z, z) is increasing along hyperbolic radii in D which 
intersect the horizontal axis. That is, let f : IB ^ Q be a conformal map 
of the unit disk IB C IR^ onto the domain Q, for which /(— 1,1) = Tj is 
the axis of symmetry of Q. Then for all all t > 0, p^{t, f{zi), f{zi)) < 
/(-^s), /(•22)), where z\ = rie*^, Z2 = r2e*^, < 6 < tt and < ri < 
ra < 1. 

For a hot-spots version, which inspired Conjecture 5.1, we refer the reader 
to [13] (Theorem 1.1) and [3] (Theorem 3.1). 

As mentioned in the introduction, for the Dirichlet heat kernel in IB , the 
opposite actually holds. That is, we have the following 

Proposition 5.2. Let IB be the unit ball in IR^, d> 2, and let p^(t,x,y) be 
the heat kernel for the Laplacian in IB with Dirichlet boundary conditions. 
(Equivalently, p-^{t,x,y) are the transitions probabilities for the Brownian 
motion in IB killed on its boundary of the ball.) Fix t > 0. The (radial) 
function p^(t, x, x) decreasing as \x\ increases to 1. That is, for all t>0, 

(5.1) p^{t, X2,X2) < p^{t, Xi,Xi), 

whenever < \xi\ < \x2\ < 1. 

Remark 5.3. People often mention this result and assert that "it is clearly 
obvious by symmetrization. " However, a proof does not seem to be written 
down anywhere. We should also mention here that the first attempt for a 
proof simply based on some type of symmetrization argument, or eigenfunc- 
tion expansion, seems to rapidly fail. For completeness, we give a short proof 
here based on the celebrated "log-concavity" results of H. Brascamp and E. 
Lieb [7]. 

Proof. Setting ^'(r) = p^{t,x,x) for |x| = r we see that by [7], for any 
< A < 1, 

(5.2) ^{Xro + (1 - A)ri) > ^(ro)^*(ri)(^-^\ 
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for any vq, ri G [0, 1]. If < ri < r2 < 1, we take tq = and pick A G (0, 1) 
such that (1 — A)r2 = ri. It follows from (5.2) that 

(5.3) ^(n) > ^'(0)^^'(r2)(^"^\ 

However, it also follows from the multiple integral re-arrangement inequali- 
ties of [7] that 

^(r) = pg{t, X, x) < pg(t, 0, 0) = ^(0), 

for all < |x| < 1. Substituting this into (5.3) proves (5.1) and completes 
the proof of the proposition. □ 

Remark 5.4. For the unit interval I = (0,1), the Neumann and Dirichlet 
heat kernels are given by 

oo 

(5.4) pf{t, x,x) = l + Y, 2e~"'"'* cos2(n7rx) 

n=l 

and 

oo 

(5.5) pf{t, x,x)=Y, 2e-"'"'* sin2(n7rx), 

n=l 

respectively; see [1]. (The "2" is there to normalize the eigenf unctions in 
L^.j Differentiating (5.4) with respect to x we see that 

d d 

(5.6) —pf{t,x,x) = -—pf{t,x,x). 
ox ox 

However, the same argument used in the proof of Proposition 5.2 above 

shows that 

(5.7) pf {t, X2,X2) < pf {t, Xi,Xi), 

whenever 1/2 < xi < 2:2 < 1- This together with (5.6) shows that pf (t, x, x) 
is increasing on (1/2, 1) and decreasing on (0, 1/2) with minimum at 1/2. 

Remark 5.5. It is interesting to note here that for the interval I, 

(5.8) pf{t,x,x)+pf{t,x,x) = Ct, 
where Ct = 1 + X^^i 2e~"'^'^^* does not depend on x. 



Since the "log-concavity" result of Brascamp-Lieb holds for all convex 
domains, the above argument gives the following result for more general 
convex domains. 

Proposition 5.6. Let be a bounded convex domain in M? which is sym- 
metric relative to the x-axis. For any {x, y) ^ Q, we write Pq {t, (x, y), (x, y)) 
for the diagonal of the Dirichlet heat kernel in fl. Fix x and let a{x) = 
sup{y > : (x,y) G ft}. The function p^ [t, {x,y), {x,y)) is decreasing in 
y for y G [0, a(x)]. By symmetry, Pq (t, (x,y), (x,j/)) is increasing in y for 
y G [-a(x),0]. 
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Motivated by the fact that the heat kernel for Brownian motion condi- 
tioned to remain forever in satisfies a Neumann-type boundary condition, 
Bafiuelos and Mendez-Hernandez proved in [5] an analogue for conditioned 
Brownian motion of the hot-spots result of Jerison and Nadirashvili [11]. 
Those results and the Laugesen-Morpurgo conjecture motivate the follow- 
ing 

Conjecture 5.7. Let ^i{x) be the ground state eigenf unction for the Lapla- 
cian in the unit ball IB C IR^, d > 1, with Dirichlet boundary conditions. 
The radial function 

p^{t,x,x) 

is increasing as \x\ increases to 1. That is, for all t > 0, 

whenever < \xi\ < \x2\ < 1. 

An appropriate version of this conjecture (similar to Proposition 5.6) mo- 
tivated by the results in can be formulated for more general convex 
domains. 
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